Highly anisotropic (/-factor of two-dimensional hole systems 



O 

o 
o 

(N 
>v 



oo 
(N 



-a 
c 

o 
o 



> 

oo 
oo 

^1- 

o 
o 
o 

-I— > 



I 

c 

O 
o 



X 



R. Winkler 

Institut fur Technische Physik III, Universitdt Erlangen-Niirnberg, Staudtstr. 7, D-91058 Erlangen, Germany 

S. J. Papadakis, E. P. De Poortere, and M. Shayegan 
Department of Electrical Engineering, Princeton University, Princeton, New Jersey 08544 

(May 28, 2000) 

Coupling the spin degree of freedom to the anisotropic orbital motion of two-dimensional (2D) 
hole systems gives rise to a highly anisotropic Zeeman splitting with respect to different orientations 
of an in-plane magnetic field B relative to the crystal axes. This mechanism has no analogue in the 
bulk band structure. We obtain good, qualitative agreement between theory and experimental data, 
taken in GaAs 2D hole systems grown on (113) substrates, showing the anisotropic depopulation of 
the upper spin subband as a function of in-plane B. 
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Since the early days of two-dimensional (2D) carrier 
systems in semiconductors it has been commonly as- 
sumed that the Zeeman energy splitting, AE — g* ^igB, 
with g* the effective (/-factor and fiB the Bohr magneton, 
is independent of the direction of the external magnetic 
field B 0| . Recently, however, calculations and experi- 
ments have shown that g* can have different values for 
B applied in the direction normal to the plane of the 
2D system compared to in-plane Here we report 

calculations and experimental data for 2D holes occupy- 
ing the heavy hole subband, demonstrating that, even 
for a purely in-plane B, g* can depend strongly on the 
orientation of B with respect to the crystal axes. We 
show that the mechanism giving rise to this remarkable g- 
factor anisotropy is fundamentally different from the one 
responsible for the in-plane/out of plane anisotropy of g* 
discussed in previous work the latter is basically 

a consequence of the bulk Zeeman Hamiltonian whereas 
the in-plane anisotropy we are reporting stems from the 
combined effect of the bulk g* and the anisotropic orbital 
motion of 2D hole states in semiconductor quantum wells. 

In bulk semiconductors the motion of electrons and 
holes in the presence of spin-orbit interaction gives rise 
to a g* which is significantly modified compared to the 
free particle g- factor go — 2 (Ref. Qj). The resulting 
g* of holes (and electrons) is nearly isotropic. Com- 
monly, the isotropic part of the hole g* is denoted by 
re g. The anisotropic part, q, is typically two orders of 
magnitude smaller than re and, in the present discussion, 
is neglected completely. The smallness of q is in sharp 
contrast to the orbital motion of holes for which we have 



highly anisotropic effective masses m* reflecting the spa- 
tial anisotropy of the crystal structure. In a 2D hole 
system (2DHS) we have heavy hole (HH) subbands (z 
component of angular momentum M = ±3/2) and light 
hole (LH) subbands (M = ±1/2). In the presence of 
an in-plane B, re couples the two LH states, and the HH 
states to the LH states Q . But there is no direct coupling 
between the HH states proportional to re. Therefore, the 
authors of Refs. concluded that the Zeeman split- 

ting of HH states due to an in-plane B is suppressed. 
However, this result is correct only for quantum wells 
(QW's) grown in the crystallographic high-symmetry di- 
rections [001] and [111] as for the other growth directions 
we show that a new mechanism gives rise to a large and 
highly anisotropic Zeeman splitting with respect to differ- 
ent orientations of the in-plane magnetic field B relative 
to the crystal axes. In the following we will discuss QW's 
grown in the crystallographic [mmn] direction (with m, n 
integers) . For this purpose we use the coordinate system 
shown in Fig. |l]with 9 denoting the angle between [mmn] 
and [001]. We remark that recently QW's for 2DHS's 
have often been grown in [113] direction as this yields 
particularly high hole mobilities ]gjJTo[] . 

We describe the hole subband states by means of 
the 4x4 Luttinger Hamiltonian ||. For the in-plane 
B = (B x ,B y ,0) we use the vector potential A = 
{zB y , —zB x ,0). Treating A and B by means of degener- 
ate perturbation theory we obtain in second order for g* 
at the bottom of the HH subbands in an infinitely deep 
rectangular QW 
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with 

7™ = -7i + 2[(l-a) 72 +a7 3 ] (2a) 

7z LH = -71 - 2 [(1 - a) 72 + a 73 ] (2b) 

a = sin 2 (6l)[3- |sin 2 (e»)]. (2c) 

Here 71, 72, and 73, are the Luttinger parameters ||] 
and 7™ and 7^ H are the reciprocal effective masses in 
z-direction in the axial approximation for the HH and 
LH subbands, respectively In Eq. (1) the factor 

(73 — 72) characterizes the anisotropy of the bulk valence 
band ]12[| . Hence, the anisotropic g-factor (1) is due to 
the combined effect of the isotropic bulk g-factor k and 
the anisotropy of the valence band. It disappears in the 
axial limit (73 — 72) = 0. The Zeeman splitting discussed 
here has no analogue in the bulk band structure. It is a 
unique feature of 2D systems. Therefore, Eq. (1) is fun- 
damentally different from the anisotropy of g* discussed 
in Refs. as the latter originates in properties of the 
bulk Zeeman Hamiltonian. Our sign convention for g* 
used in Eq. (1) corresponds to the dominant spinor com- 
ponent of the multicomponent eigenstates (using a basis 
of angular momentum cigenfunctions with quantization 
axis in the direction of B). Note that in unstrained QW's 
the topmost subband is the HH1 subband. 

For LH subbands in an in-plane B as well as for HH 
and LH subbands in a perpendicular B, g* contains terms 
similar to Eq. (1). However, the dominant contribution is 
given by the bulk g-factor k. For LH subbands in an in- 
plane B we have basically g^ — 4k, while for a perpen- 
dicular B we have gf H = 6k and g^ H = 2k (Refs. |§|§). 
Yet, this implies a remarkable difference |^| compared 
with Eq. (1). 

In Fig. we show the anisotropic g* of the HH1 sub- 
band for a 200 A wide GaAsZAlo.3Gao.7As QW as a func- 
tion of the angle 9. The analytical expressions (1) (dot- 
ted and dashed-dotted lines) are in very good agreement 
with the more accurate results obtained by means of a 
numerical diagonalization [Q,[l3| of the Luttinger Hamil- 
tonian (solid and dashed lines). Figure |^ demonstrates 
that g* can be very anisotropic |14[| . For example, for 
the growth direction [113], g* is about a factor of 4 larger 



when B || [332] compared to when B || [110]. Moreover, 
the sign of g HH - , is opposite to the sign of oS^, . 

[nn{2m)\ L lltU J 

Equation (1) is applicable to a wide range of cubic 
semiconductors with results qualitatively very similar to 
Fig. |^. In particular, the relative anisotropy 

g [Tio] _ sin(<?) /„s 

HH - I — W 

y [nn(2m)] W 4 - 3 Sm (9) 

is independent of the material-specific parameters 7 j and 
k. This remarkable result can be traced back to the fact 
that the anisotropy for different directions 9 in k-space 
is always characterized by the single parameter (73 — 72) 
(Ref. [Q). Note that for QW's based on narrow-gap 
semiconductors we have a larger k and smaller effective 
masses so that for these materials the absolute values 
of g* are significantly larger than g* of GaAs shown in 
Fig. ||, but the g* anisotropy is still given by Eq. (||) and 
depends only on 9. The Zeeman splitting can be even fur- 
ther enhanced if one uses semimagnetic semiconductors 
containing, e.g., Mn. For these materials the structure 
of the Hamiltonian is identical to the conventional Lut- 
tinger Hamiltonian in the presence of a magnetic field 
with k replaced by the effective <?-factor due to the para- 
magnetic exchange interaction |^||. Therefore Eq. (1) 
and the g* anisotropy [Eq. (^)] are readily applicable to 
semimagnetic materials, also. 

In a Taylor expansion of the Zeeman splitting, AE(B), 
g* (times /is) is the prefactor for the lowest order term 
linear in B. Often terms of higher order in B are ne- 
glected because of their relevant insignificance. An in- 
teresting feature of Fig. |^ is the vanishing of g* for the 
high-symmetry growth directions [001] and [111] (Ref. 

For the 2DHS discussed in Fig. || this results in a 
splitting AE which at B = 1 T is more than 2 orders of 
magnitude smaller than AE for growth directions [113] 
and [110]. For these high-symmetry directions AE is pro- 
portional to B 3 (Ref. |15| ) . In second order perturbation 
theory we obtain for the HH1 subband of an infinitely 
deep rectangular QW of width d grown in [001] direction 



A£=Mf(^) V / 72 2 cos2(2^+ 7 3 2 Bin 2 (20) 



Here cj> is the angle between the in-plane B and the [100] 
axis. The first term in the square bracket stems from kp 
coupling between the subbands HH1 and LH1, and the 
second term is due to coupling between HH1 and LH3. 
We get similar, though somewhat longer expressions for 
growth direction [111]. It is remarkable that we have a 
nonzero AE even in the limit k = 0. This can be un- 
derstood as follows: The 4x4 Luttinger Hamiltonian ||] 
which is underlying our calculations corresponds to an 



4k (tt 2 - 6) 27 7 3 

^HH ^LH ^HH q^/LH 



infinitely large spin-orbit splitting between the topmost 
valence band Tg and the split-off band Ty. Therefore 
spin-orbit interaction is not explicitely visible in our re- 
sults, though, similar to Zeeman splitting in bulk semi- 
conductors 0, Eqs. (1) and (||) are a consequence of 
spin-orbit interaction. In 2D systems the motion of elec- 
trons and holes in the presence of this interaction can 
give rise to a Zeeman splitting even without a bulk g* . 
We remark that in a parabolic QW the in-plane g* of the 
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HH subbands also contains such terms independent of k. 
We have here a 2D analogue of Roth's famous formula J7 ] 
for the electron bulk g*. Finally we note that, unlike |fl] 
Eq. (1), AE in Eq. (^) increases proportional to d 4 , i.e., 
Zceman splitting is most efficiently suppressed in narrow 
QW's. 

Now we will show that the anisotropy of g* can be 
probed experimentally by measuring the magnetoresis- 
tance of a high-mobility 2DHS as a function of in-plane 
B. The samples are 200 A wide Si- modulation doped 
GaAs QW's grown on (113)A GaAs substrates. These 
samples exhibit a mobility anisotropy believed to be due 
to an anisotropic surface morphology |]l6|] . They are pat- 
terned with an L-shaped Hall bar to allow simultaneous 
measurements of the resistivity along the [332] and [110] 
directions. Front and back gates are used to control the 
2D density in the QW and the perpendicular electric field 
which characterizes the asymmetry of the sample [juJ . 

The left two panels of Fig. || show the resistivity p 
measured as a function of in-plane B for three differ- 
ent densities and different relative orientations of B and 
current /. For easier comparison we have plotted the 
fractional change p(B)/p{B = 0). Apart from an over- 
all positive magnetoresistance these curves show a broad 
feature consisting of an inflection point followed by a re- 
duction in slope followed by another inflection point. In 
Fig. U we have placed arrows between the two inflection 
points at a value of B we call B* . Similar, though sharper 
features have been observed in systems with several oc- 
cupied confinement subbands when a subband is depop- 
ulated by means of an in- plane B Q. Wc propose that 
the magnetoresistance feature at B* in Fig. |^ is related to 
a spin-subband depopulation and the resulting changes in 
subband mobility and intersubband scattering as the in- 
plane B is increased. Note that in each panel of Fig. ||, B* 
is the same for both current directions even though the 
magnetoresistance is very different. This implies that B* 
depends on parameters which do not depend on current 
direction. This supports our hypothesis, as spin-subband 
depopulation should not depend on the direction of cur- 
rent in the sample. 

Our interpretation of B* is obviously consistent with 
B* in Fig. [3] becoming larger with increasing density. It 
is remarkable that B* for the B || [332] traces is about 
4 T smaller than for the [11 0] traces, regardless of the 
/ direction. We associate this with the anisotropy of 
the in-plane g*. This interpretation is validated by our 
self-consistently calculated J|,[ll],[l3| results for the den- 
sity N+ of the upper spin subband as a function of B 1 
shown in the right panel of Fig. |^. The density N + de- 
creases much faster for B || [332] than for B || [110], in 
agreement with Fig. |^. We have further support for our 
interpretation of B* from experiments where we increase 
the asymmetry of the confining potential by means of 
the front and back gates while keeping the 2D density 
fixed. We observe an increase in B*, in agreement with 
the results of our self-consistent calculations. 



One might ask whether the data in Fig. [3| actually 
could be summarized by a single value of g* for each 
trace. Unfortunately, this is not possible because, due to 
the complicated band structure of holes, g* depends on 
energy E, and we are averaging over g*(E) for E between 
the subband edge and the Fermi energy. The importance 
of this effect can be readily deduced from the right panel 
of Fig. ||, as we would have straight lines for N + (B) if g* 
(and the effective mass m*) were not dependent on E. 

In Fig. H the measured B* is significantly smaller than 
the calculated B for a complete depopulation of the upper 
spin-subband. We note that for our low-density samples 
it can be expected that g* is enhanced due to the ex- 
change interaction jl]Jl8|| . This effect was not taken into 
account in our self-consistent calculations. The qualita- 
tive agreement between the experimental data and the 
calculations, however, implies that these many-particle 
effects do not affect the anisotropy of g*. 

The large anisotropy of the Zeeman splitting in 2DHS's 
offers many possible device applications. In a poly crys- 
talline material, e.g., one could alter the degree of spin 
polarization in different domains by changing the direc- 
tion of the external B. As can be seen in Fig. ||, because 
of the sign reversal of g* it is even possible to have differ- 
ent domains with opposite spin polarization for a given 
direction of B. Recently, there has been a growing in- 
terest in controlling the spin degree of freedom for quan- 
tum computing and spin electronics. In Ref. JlSj ] the au- 
thors have sketched a quantum device which makes use 
of the spatial variation of g* in layered semiconductor 
structures made of, e.g., Al a Jn !; Gai_. ! ;_ !/ As. However, 
the authors have estimated that a substantial change in 
g* requires a fairly large electric field of the order of 
100 kV/cm. Oestreich et al. |2(| and Fiederling et al. 
plf have suggested a spin aligner based on semimagnetic 
semiconductors. Here the (/-factor anisotropy of 2DHS's 
provides a powerful additional degree of freedom for en- 
gineering such devices. 

In conclusion, we have shown that coupling the spin 
degree of freedom to the anisotropic orbital motion of 
2D hole systems gives rise to a highly anisotropic Zee- 
man splitting with respect to different orientations of an 
in-plane magnetic field relative to the crystal axes. 
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FIG. 1. Coordinate system for QW's grown in [mmn] di- 
rection (z direction). Here 8 is the angle between [mmn] and 
[001], i.e., we have 8 — arccos(n/\/2m 2 -f n 2 ). The axes for 
the in-plane motion are [nn(2m)] (x) and [TlO] (y). 
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FIG. 2. Anisotropic effective g-factor g* of the HH1 sub- 
band for a 200 A wide GaAs/Alo.3Gao.7As QW as a function 
of 8, the angle between [001] and the growth direction. Re- 
sults are shown for the in-plane B along the [nn(2m)] and 
[TlO] directions. The solid and dashed lines were obtained by 
means of a numerical diagonalization of the Luttinger Hamil- 
tonian. The dotted and dashed-dotted lines were obtained by 
means of Eq. (1). 
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FIG. 3. Left and central panels: Fractional change in resis- 
tivity p(B)/p(0) due to an in-plane B, measured at T = 0.3 K 
in a GaAs 2D hole system grown on a (113) substrate, for 
different 2D densities as indicated. The arrows mark B* as 
defined in the text. Right panel: Calculated density A+ in 
the upper spin subband as a function of B. 
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